Search for: All records

Abstract RNA dependent RNA polymerase (RdRp), is an essential in the RNA replication within the life cycle of the severely acute respiratory coronavirus-2 (SARS-CoV-2), causing the deadly respiratory induced sickness COVID-19. Remdesivir is a prodrug that has seen some success in inhibiting this enzyme, however there is still the pressing need for effective alternatives. In this study, we present the discovery of four non-nucleoside small molecules that bind favorably to SARS-CoV-2 RdRp over the active form of the popular drug remdesivir (RTP) and adenosine triphosphate (ATP) by utilizing high-throughput virtual screening (HTVS) against the vast ZINC compound database coupled with extensive molecular dynamics (MD) simulations. After post-trajectory analysis, we found that the simulations of complexes containing both ATP and RTP remained stable for the duration of their trajectories. Additionally, it was revealed that the phosphate tail of RTP was stabilized by both the positive amino acid pocket and magnesium ions near the entry channel of RdRp which includes residues K551, R553, R555 and K621. It was also found that residues D623, D760, and N691 further stabilized the ribose portion of RTP with U10 on the template RNA strand forming hydrogen pairs with the adenosine motif. Using these models of RdRp, we employed them to screen the ZINC database of ~ 17 million molecules. Using docking and drug properties scoring, we narrowed down our selection to fourteen candidates. These were subjected to 200 ns simulations each underwent free energy calculations. We identified four hit compounds from the ZINC database that have similar binding poses to RTP while possessing lower overall binding free energies, with ZINC097971592 having a binding free energy two times lower than RTP. 

Hepatitis C Virus (HCV) is reported to develop GS-461203 resistance because of multiple mutations within the RNA-dependent RNA Polymerase (RdRp) of HCV. The lack of a high-resolution structure of these RdRp mutants in complex with GS-461203 hinders efforts to understand the drug resistance. Here we decipher the binding differences of GS-461203 in the wild type and mutated systems T179A or M289L of HCV RdRp Genotype 2a using homology modeling, molecular docking, and molecular dynamics simulation. Key residues responsible for GS-461203 binding were identified to be Arg48, Arg158, Asp318, Asp319, and Asp220, and that mutations T179A or M289L have caused conformational changes of GS-461203 in the RdRp active site. The affinities of GS-461203 were reduced in T179A system, but it became slightly stronger in the M289L system. Furthermore, we designed two new analogues of GS-461203 which encouragingly induced more stable interactions than GS-461203, and thus resulted in much better binding energies. This present study reveals how a single mutation, T179A or M289L, will modulate GS-461203 binding in HCV RdRp Genotype 2a, while introducing two novel analogues to overcome the drug resistance which may be good candidate for further experimental verification. 

The novel coronavirus disease 19 (Covid-19) which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a pandemic across the world, which necessitate the need for the antiviral drug discovery. One of the potential protein targets for coronavirus treatment is RNA-dependent RNA polymerase. It is the key enzyme in the viral replication machinery, and it does not exist in human beings, therefore its targeting has been considered as a strategic approach. Here we describe the identification of potential hits from Indonesian Herbal and ZINC databases. The pharmacophore modeling was employed followed by molecular docking and dynamics simulation for 40 ns. 151 and 14480 hit molecules were retrieved from Indonesian herbal and ZINC databases, respectively. Three hits that were selected based on the structural analysis were stable during 40 ns, while binding energy prediction further implied that ZINC1529045114, ZINC169730811, and 9-Ribosyl-trans-zeatin had tighter binding affinities compared to Remdesivir. The ZINC169730811 had the strongest affinity toward RdRp compared to the other two hits including Remdesivir and its binding was corroborated by electrostatic, van der Waals, and nonpolar contribution for solvation energies. The present study offers three hits showing tighter binding to RdRp based on MM-PBSA binding energy prediction for further experimental verification.